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We consider the Schrodinger functional with staggered one-component fermions on a fine lattice 
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component spinors, two different set-ups are proposed, corresponding to the coarse lattice having 
size (L/2fl)^ X (T' /2a), with T' = T ±a. The continuum limit is then defined at fixed T' /L. 
Both cases have previously been investigated in the pure gauge theory. Here we define fermionic 
correlation functions and study their approach to the continuum limit at tree-level of perturbation 
theory. 
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1. Introduction 

Establishing a quantitative connection between the low energy and the perturbative regimes of 
QCD is one of the primary tasks for any attempt to solve QCD quantitatively. For lattice QCD, the 
main problem consists in the large scale differences involved, which cannot be resolved on a single 
finite lattice. A solution to this "non-perturbative renormalisation problem" has been proposed 
a while ago [|I|] and amounts to apply recursive finite size scaling techniques to the renormalised 
parameters and operators in a suitable renormalisation scheme. The Schrodinger functional (SF) [§, 
|3|] gives rise to a class of such schemes with a number of technical advantages. In QCD with zero 
and two quark flavours, the running coupling ||^, 0, ^, 0] and quark masses[||, ||, ^], as well as 
a range of composite operators have been studied. Note that the final results are obtained in the 
continuum limit and thus independent of the details of the lattice regularisation. Most results have 
been obtained using the implementation of the SF in QCD with Wilson type quarks [||]. However, 
for applications to QCD with four quark flavours, or for QCD like theories with multiples of four 
fermion flavours staggered fermions appear to be a natural alternative. Interesting universality tests 
could be devised and one may expect a better control over the continuum limit. Hence, in view of 
applications to four-flavour QCD we here revisit the implementation of the Schrodinger functional 
for staggered quarks, which has previously been studied in [[I^, Its applications have so far 
been limited to studies of the running coupling for QCD-like theories with eight, twelve and sixteen 
fermion flavours [11, 12 1. As noticed in [[l0|, pi]], the time extent of the lattice, T /a, needs to be 



odd with staggered quarks, whereas the spatial lattice directions must have even extent, L/a. This 
makes it impossible to set L = T exactly, and one needs to deal with the resulting 0(a) effects. In 



order to cancel those. Heller [[11|] proposed to average results for the gauge coupling obtained in 
two separate simulations with T =L±a. While this seems to work out for the SF coupling, at least 
to one-loop order in perturbation theory, it is less clear how to proceed in the case of fermionic 
correlation functions. In particular, one needs to discuss how to reconstruct the four-component 
spinors in both cases, only one of which was considered in |[TT|]. Ideally, one would like to avoid 



the averaging procedure altogether, and it has been shown in the pure gauge theory how this can be 



achieved by redefining the approach to the continuum limit [|13[]. 

This writeup is organised as follows. We start by reviewing the basics of the Schrodinger 
Functional and the definition of fermionic correlation functions. Next we reconstruct the action in 
terms of the four-component spinors for both cases, T' = T±a. A chiral rotation is then carried out 
to recover the standard Schrodinger Functional boundary conditions for the fermionic fields. We 
show the results of the computation of the correlation functions at tree level of perturbation theory 
and we finish with an outlook to future work. 

2. The Schrodinger Functional and correlation functions 

The Schrodinger Functional is the Euclidean path integral of QCD on a hyper cylinder as 
space-time manifold. Dirichlet boundary conditions are imposed at Euclidean times xq = 0, T, 
while all fields are L-periodic in the spatial directions. For the fermionic fields one sets 
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where P± = ^(1 ± 7o)- Using a continuum notation, the spatial gauge field components satisfy the 
conditions, 



The Schrodinger Functional can then be regarded as a functional of the boundary fields. 



-S[A,v/,V/] 



and expectation values of any product of fields ^, are defined by, 

1 



(2.2) 



(2.3) 



(2.4) 



p=p'=0;p=p'=0 



Note that observables may contain quark and antiquark fields at the boundaries, by including deriva- 
tives with respect to the fermionic boundary fields, (^(y) = gp^(y) i C (y) = ^^'^ analogously 
for (^'(y), (^'(y). Provided the gauge boundary fields are taken to be spatially constant, one may 
obtain gauge invariant quark bilinear sources at the boundaries, such as 



(2.5) 



where x" is a flavour matrix and both the quark and anti-quark fields are projected to zero momen- 
tum. Using such sources the simplest fermionic correlation functions for the axial vector current 



and density take the form | |14| ] 

ft{y,) = -{Al{y)&'), ff{yo) = -{P^iy)^'), ff = (2.6) 
Note that with an exact flavour symmetry all correlation functions would be proportional to 5"^. 

3. Reconstruction of the four component spinors 

3.1 Case r = r -a 

This case is the one already discussed in[|To|, pT] |. The four-component spinors reside in a 
coarse lattice with lattice spacing a = 2a. In Figure |l], the thin lines correspond to the fine lattice, 
and the dots represent the points of the coarse lattice where the reconstructed fermions live. The 
variable y refers to the points in the coarse lattice, and x to the fine lattice, and they are related 
hy X = 2y + a^, with taking values in {0,1}. Introducing the transition fields X£,{y) = X{x)> 

Dirichlet 

Q+V'T — — \ — T — I Q-^ 



^0 1 2 

a 

Figure 1: Reconstruction of the spinors onaT = L + a lattice. 
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{y) = X{^)- the transformation is specified as 



(3.1) 



with = \y^y\y^y^ ■ in Figure |l|, the one-component fermionic fields which constitute a re- 
constructed quark field are the ones contained in the circles. Note that the Dirichlet boundary 
conditions at = 0, T , imply a projection onto half of the components of the reconstructed quark 
field. Labelling the (hermitean) flavour matrices by their 7-matrix structure, e.g. = 7^,1^5 = 
i{y^y5Y and denoting the symmetric derivative by and the second derivative by A^, the 
boundary conditions read 



e+VA(0,y) = p(y), Q-YiT',y) = p'(y), 
VA(0,y)G+ = p(y), Y{T',y)Q- = p'(y), 



(3.2) 



with projectors Q± = ^(1 it /7075T05). For homogeneous boundary conditions, and with all fields 
at times xq < and xq > T' set to zero, the reconstructed action takes the form. 



.(,v=-l) 



yo=o y M 



yfid^i+i-ysT^s^fi 



viy)- 



(3.3) 



The usual SF boundary conditions can be recovered by performing a chiral rotation of the fermionic 
fields, 

\l/{y)=R{a)Y{y), Y'{y) = ¥{y)R{oc), /?(«)= exp(/f 75105)- (3.4) 

For « = I the boundary conditions become the usual ones (p]), due to R{j)Q±R^^{j) = P±. For 
homogeneous boundary conditions, the action in the standard SF basis takes the form. 



.(.v=-l) 



T' 



y=0 y 



Y,yk^k + yQdo + -Ao 
k ^ 



Viy), 



(3.5) 



with &k=~dk + i^ykys'^ksAk- 
3.2 Case T' = T + a 

Here, we distinguish two alternative ways of reconstructing the fermions, as illustrated in 
Figure ^. We have labelled the two reconstructions with 5 = 1*, according to the sign in front of 
^0 in Eq. ( |3^ . 



s = l+ 

xo = 2yo-a + a^o, x = 2y + a^, 

¥aa{y) = 5£(f^)a«Z^(3'), 



3lZ^(3')(fp.,a, 



Vaaiy) -- 

f^ = i(-l)^°7?,7f'7?7? 



s = r 

xo = 2yQ-a^o, x = 2y + a^, 

Vaaiy) = \J^i^^)aaX^iy), 
Vaaiy) = \Y.X^iy)i'[^l)aa 



(3.6) 
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Figure 2: Reconstruction of the spinors on a T = L — a lattice. Left s ^ 1 +, right, 5=1" 



The interpretation of Figure |^ is the same as Figure [IJ Eqs. (3.6) specify how to reconstruct the 
four-component fermions in both cases. Depending on the reconstruction, the boundary conditions 
are different, and therefore the chiral rotations needed to restore the usual SF b.c.'s differ, too: 

5 = 1+ * = 1" 

G-i/^(0,y)=p, Q+Y{0,y)=p', Q+¥{0,y)=p, Q-Y{0,y)=p', 

¥iO,y)Q-=P, w'{0,y)Q+=p, i^(0,y)e+ = P v/(0,y)e_=p, (3.7) 

However, once rotated to the standard SF basis, the action for both cases s = 1+ is the same. 



yo,y 



k ^ 



(3.8) 



3.3 Symmetries of the SF with staggered fermions 



The symmetries of the SF with staggered quarks have been summarised by Heller []1 IQ. We 
here just wish to emphasise that the flavour and chiral symmetries refer to a particular basis. As we 
have seen, the boundary conditions may depend on the way the four-spinors are reconstructed. It is 
only after performing a chiral non-singlet rotation that the standard SF is recovered. In this basis, 
the usual axial U(l) symmetry of staggered quarks becomes a flavour symmetry. More precisely, 
the transformation 



¥iy) 



¥iy), 



¥iy) ¥iy)e 



when rotated into the SF basis, becomes a continuous flavour symmetry with generator To> 

^(3;) ^ e'l^'Wiy), Y'iy) - v'iy)e-'^'°- 

Furthermore we note that spatial translations by a unit a on the fine lattice. 

¥{y) T:kY{y) + mQ+ik)dkY{y) ¥{y) ¥{y)'^k + a¥{y)dk'Ckd-^^ 



+ ' 



(3.9) 



(3.10) 



(3.11) 



with Q)^' = ^(1 lb ykJs'Cks)', correspond to a discrete subgroup of flavour symmetry in the continuum 
limit. 

4. Correlation functions at tree level 

To evaluate the correlation functions we first integrate over the quark fields. The expectation 
value assumes the form (1^) = )g> where ()g denotes the gauge field average. We have deter- 
mined the free quark propagator both analytically and numerically and may therefore compute the 
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correlation functions to tree level. The remaining chiral symmetry of Eq. (3.9) becoming a flavour 
symmetry with generator Tq in the standard SF basis, disconnected diagrams in the computation of 
fA,fp,fi are forbidden if we choose those flavour matrices which anticommute with Tq. For these 
matrices, the correlation functions /^^(jo) reads, 

ffiyo) = a' lj(tr ([C(y'')V/^(3')]F7o75Tl'/WC(y')]F75^'))^, (4.1) 

and analogous expressions are obtained in the other cases. The continuum values of at tree 
level with vanishing background field take the form, d"^fx, with 

^^^"^'/'^^-cosh^Tvse)' ^^^"^'/'^^cosh^Vse)' ^' = cosJiV3ey ^^"'^ 

where is a phase factor coming from the generalised boundary conditions, i.e. + Lk) = 
e'^\ir{y):¥iy + ^^) = Viy)^ '^- including the correct tree level boundary counterterm, the results 
obtained are accurate up to O(a^) for /p,/i and 0(a) for /a. 

5. Fermionic 0(a) improvement 
5.1 Infinite volume 

Close to the continuum limit, the lattice theory may be described in terms of a local effective 
theory with action [15], 

5eff = So + aSi + a^S2 + ..., Sk = J d^y^k {y) (5.1) 

The apparent 0{a) contributions on the infinite lattice are fixed by the shift symmetry, since it is the 
combined expression ^^Y^ that is invariant under this transformation and not the usual kinetic term 
alone As was pointed out in |16], there are no invariant dimension 5 operators, so no counterterms 
can be added. The standard procedure to eliminate the apparent 0(a) terms consists in defining 
improved field, 

w\y) = ¥iy) + lUQ+^-Q^-ydMy), 



V 



v'iy) = m + lLv^viym?-Q^-^), (5.2) 



4 ^ 

in terms of which one finds, 

SsQ = a'l^Y'{y)Y^~d^Y'{y) + 0{a^). (5.3) 

5.2 0{a) effects from the boundaries 

In the SF framework, additional renormalisations and 0(a) cutoff effects may arise from the 
very presence of the boundaries. Taking the symmetries into account, we arrive at the conclusion 
that there is only one possible dimension 3 operator, y'y'- This is the same as encountered for 
Wilson quarks, and can thus be absorbed in a multiplicative renormalisation of the quark and anti- 
quark fields at the boundaries. 
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In the case of dimension 4 operators, we obtain again the same result as for Wilson quarks. 
However, when using the equations of motion, we here prefer a different choice for the counterterm 
action, namely 

55f [f/, VA, VA] = £ { (Ji - 1 ) [4,1 + + (^2 - 1 ) [4,2 + ^,,2] } , (5.4) 
y 

4,1 = ¥'{o,y)P+rk^^kv'{o,y), 4.2 = Piy)rk2^kp{y), (5.5) 
4,1 = v{T,y)P-rk^k¥{T,y), 4,2 = p'(y)u^^p'(y), (5.6) 

The coefficients di^2 have a perturbation expansion in powers of g^. We have determined di at tree 
level, 

, =l±j. (5.7) 

T'=T±a 4 

6. Conclusions 



We have reconstructed the four-component spinors in the SF with staggered quarks, for both 

cases T' = T ±a, computed the free propagator and the correlation functions /a,/p,/i at tree 
level. The implementation of 0(a) improvement is work in progress. Once it is fully understood, 
we plan to trace the running of the SF coupling and the quark mass in four-flavour QCD 
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